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I.  Introduction 
A.  Overview 

The  purpose  of  this  report  is  to  discuss  solution  methods  for  several 
problems  arising  in  inventory  provisioning.  Two  related  types  of  problems  are 
considered: 

1)  optimization  of  a  performance  measure  subject  to  a  budget 
constraint 

2)  minimization  of  cost  subject  to  a  constraint  on  performance. 

The  performance  measures  considered  are  Supply  Material  Availability  (SMA), 
Mean  Supply  Response  Time  (MSRT),  and  "Pseudo-Availability"  (PA).  The  basic 
approach  to  these  problems  is  dynamic  programming. 

First  the  problems  are  fornulated  and  solved  using  one  recursive 
technique.  Then  a  more  efficient  recursion  is  presented  and  discussed  in  the 
context  of  maximizing  MSRT  subject  to  a  budget  constraint. 

These  same  budget  constrained  problems  are  formulated,  discussed  and 
solved  using  a  marginal  analysis  approach  in  [2]  and  [3],  That  method, 
although  fast,  does  not  guarantee  that  optimal  solutions  are  obtained.  This 
report  provides  a  method  of  obtaining  optimal  solutions  and  thus  provides  a 
means  of  evaluating  heuristic  methods.  In  addition,  this  report  shows  how  the 
performance  constrained  budget  minimization  problems  can  be  solved  directly. 

By  contrast,  reference  [2]  addresses  this  problem  using  generalized  Lagrange 
Multipliers  or  by  solving  the  budget  constrained  performance  problem 
repeatedly  for  various  budget  levels. 

The  second  recursion  presented  here  Is  considerably  faster  than  the  first 
and  It  also  guarantees  an  optimal  solution  in  some  cases.  It  Is  competitive 
In  speed  with  the  marginal  analysis  method  for  small  and  medium-sized 


problems,  but  it  is  inefficient  for  problems  with  a  large  number  of  items  or 
large  budget  values.  Its  virtue  lies  in  its  ability  to  get  exact  solutions  to 
medium-sized  problems. 

B.  Problem  Formulations,  budget  constrained 

This  report  considers  three  specific  budget  constrained  optimization 
problems  arising  in  Inventory  provisioning.  These  problems  are  formulated  and 
discussed  in  detail  elsewhere  [2]  and  will  be  stated  here  without  extensive 
explanation. 

Following  [2]  we  let 

n  *  the  total  number  of  items  considered  for  provisioning 
Cj  *  the  unit  cost  of  item  1 

*  the  essentiality  code  for  item  type  i 

*  the  demand  rate  for  item  i 

Tj  *  the  procurement  leadtime  for  item  i 

=  the  number  of  items  of  type  i  provided  (the  decision 
variables) 

Zi(Si)  *  the  performance  measure  for  item  i  when  Si  units  are 
stocked 

*  zi<si-i>  -  zi<si> 

Pj(Xj)  3  probability  that  demand  for  item  i  Is  x^  during  the 
provisioning  Interval 

P^xj)  3  cumulative  probability  of  x^  or  fewer  demands  during 
a  provisioning  Interval 

MTTR^  3  mean  time  to  repair  or  replace  item  1 
MTBF^  *  mean  time  between  failures  3  1/A^ 

MSRTj(Sj)  3  mean  supply  response  time  when  units  are  stocked. 
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The  three  budget  constrained  problems  considered  are: 

al)  maximize  Supply  Material  Availability  (SMA),  defined  as 


SHA(Sl . S„)  •  E,  1,  T,  Z^IS,)/^  E,  X,  T, 

Z,(l)(St)  -  (1  -  P,(S,))  ♦  (S,  -  X,  rt)(l  -  Pt(S1))/X1  T 
bl)  minimize  Mean  Supply  Response  Time  (MSRT),  defined  as 

MSRT(Slt...,Sn)  *  ^  X^  T^  Z-  ^Si  )/^j  Ei  Ai  T1 

12)  W1*, 

z,‘2)<V  .  a  -  WHyz  -  St/x,  ♦— y— -> 

^  i  i 

♦  P1(S1)(X1  T1  -  Si)/2Xi  . 

cl)  maximize  Pseudo-Availability  (PA),  defined  as 

PA(S,,....S  )  -  I  zP^SJ 
1  n  w  i  i 

zj3)(S  )  •  MTBF1/(MTBF1  +  MTTR^  +  z}2^))  . 
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II.  Solution  Method  and  Examples 


A.  Dynamic  Programming  Approach 

The  computer  program  used  to  solve  these  problems  Is  DP4,  a  general 
purpose  program  for  performing  dynamic  programming  tabular  computations.  Here 
we  will  describe  the  general  nature  of  that  program  and  those  elements 
required  to  tailor  it  for  use  in  the  problems  considered  in  this  report. 

The  DP4  program  deals  with  a  problem  consisting  of  n  related  stages 
each  of  which  is  characterized  as  shown  in  figure  1. 


*1  - *  '  — •  xi.i  •  MvV 

^ “ 
ri(xi  •‘M 

Figure  1.  A  single-stage  decision  problem. 

In  figure  1 

Xj  is  the  "state"  variable 

dj  is  the  decision  variable 

r^  is  the  stage  return  function 

tj  is  the  stage  transformation  function. 

In  the  overall  problem  consisting  of  n  stages  the  output  state  variable  from 
stage  1  ,  namely  X^_j  ,  is  the  input  to  stage  1-1  .  Thus,  the  n  stage 
problem  can  be  pictured  as  in  figure  2. 


Figure  2.  n-stage  decision  problem. 


The  state  variable  can  easily  be  understood  in  the  context  of 
the  budget  constrained  problems  as  the  amount  of  resource  (money)  remaining  to 
be  allocated  to  stages  i,  1-1, ...,1  .  The  stages,  of  course,  correspond  to 
the  items  in  the  inventory  problems. 

At  each  stage  i  a  decision  d^  must  be  made.  The  decision  has  two 
effects.  First,  it  yields  a  return  r^  ,  the  performance  measure  associated 
with  the  current  item.  Second,  it  yields  a  value  of  X^_1  which  serves  as 
the  input  to  the  remainder  of  the  decision  process.  The  decision  d.  must  be 
made  with  consideration  both  for  the  immediate  return  r.  and  the  future 
state  Xj_j  .  The  overall  problem  is  to  make  the  series  of  decisions 
dn,...,d^  to  optimize  some  function  of  the  individual  stage  returns. 

In  the  problem  (al),  where  is  the  decision  variable,  we  can  let 
the  return  functions  be 


1(Xi,  St)  -  Ei  Xi  T.  Z1(1)(Si  )/iI1E1  A1  T.  i  -  1 . n 


and  the  stage  transformation  functions  be 


Xn  *  8 

X1-l  *  W  V  *  X1  "  C1  S1  1  *  1 . n  • 

The  overall  return  function  Is  the  sum  of  the  individual  return 
functions.  Namely, 


SMA(Sj . Sn)  -  ^  r^,  . 

The  object  remains  to  select  Sj,...,Sn  to  optimize  this  return. 
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We  let  f^X^)  3  the  optimal  total  return  from  stages  i,  i-l,...,l 

given  that  we  enter  stage  i  with  state  variable  X^. 
Then  we  can  write  the  recursive  equations  for  this  optimization  as 


max{r . 
S1 

(X.,  S.)  ♦fi.1(X1_l 

s.t. 

Xi-1  "  Xi  “  Ci  Si 

and 

0  <  Si  <  x./c. 

and 

Si  *  integer 

for  i  3  2,... ,n  . 

The  equation  at  stage  1  is 


fj{Xj)  s  roex  ^(Xj,  Sj) 
S1 


s.t.  0<Sl<Kl/Cl 


and  S1  integer 


The  program  DP4  performs  this  optimization  provided  the  user  supplies 
the  following  subroutines  and  data. 

Required  subroutines 

1.  STGRET  -  this  subroutine  defines  the  function  r^(Xj,  ) 

2.  TRANFM  -  defines  the  stage  transformation  function  t^(X^,  d^) 

3.  DLIMIT  -  defines  the  range  of  decision  values  d,.  which  can  be 

considered  for  the  particular  value  of  X^  under 
consideration 

4.  STORE  -  allows  the  Input  of  constants  to  be  used  In  the  other 

subroutines. 
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Required  Data 

1.  n  -  the  number  of  stages 

2.  For  each  stage  i 

XLOW  -  the  lowest  value  of  X..  to  consider 

XHI6H  -  the  highest  value  of  X..  to  consider 

DELX  -  the  increment  for  X,. 

XMODE  -  tells  whether  to  maximize  or  minimize 

XSTAGE  -  tells  how  this  stage  return  relates  to  lower  numoered  stage 
returns  (sum,  product). 

The  methodology  is  essentially  the  same  for  the  performance  con¬ 
strained  problems.  There  the  return  functions  d-j )  =  c.d^  .  The 

state  variable  X.  is  interpreted  as  the  portion  of  the  performance  measure 
to  be  attributed  to  stages  .  The  stage  transformation  functions  in 

problem  (a2)  and  (b2)  are 

X.  ,  =  X.  -  Z.(S.)  . 

i-l  i  l  '  i  ‘ 

In  problem  (c2)  the  stage  transformation  is 

xi-i  *  Vz(3)(V  • 

The  program  DP4  and  the  subroutines  are  shown  in  Appendix  A.  The 
subroutines  are  written  to  solve  any  of  the  problems  (al),  (bl),  (cl)  or  (a2), 
(b2),  (c2).  Thus  they  involve  conpllcations  not  needed  for  solving  just  one 
of  these  problems. 

B.  Examples 

Several  example  problems  were  solved  to  Illustrate  the  approach  discussed 
here.  All  of  the  problems  Involved  n  ■  10  Items  and  all  used  the  data  shown 
in  table  1. 
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Data 


1.  The  budget 
Tables  2, 
(bl),  and  (cl) 


n 

X 

Time 

Cost 

MTTR 

E 

1 

5.0 

1.0 

1.0 

.0137 

1.0 

2 

2.0 

1.0 

2.0 

.0274 

1.0 

3 

3.0 

1.0 

5.0 

.0137 

1.0 

4 

5.0 

1.0 

10.0 

.0822 

1.0 

5 

10.0 

1.0 

20.0 

.0274 

1.0 

6 

25.0 

1.0 

5.0 

.0027 

1.0 

7 

1.0 

1.0 

1.0 

.0054 

1.0 

8 

1.0 

1.0 

100.0 

.0411 

3.0 

9 

0.5 

1.0 

50.0 

.0082 

1.0 

10 

2.0 

1.0 

10.0 

.1370 

3.0 

Table  1:  Data  for  Examples 
constrained  problems. 

,  and  4  sumnarize  the  solutions  for  the  example  problems  (al). 
These  are  the  budget  constrained  problems. 
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B  »  Budget  * 

max  SMA  s 

300 

.750654 

295 

. 746466 

290 

.741160 

285 

.736563 

280 

.732964 

decision  a 

4 

4 

4 

4 

4. 

S2* 

4 

4 

4 

4 

4 

S3* 

4 

4 

4 

4 

4 

S4  * 

5 

5 

5 

4 

5 

S5  * 

2 

2 

2 

2 

1 

S6* 

29 

28 

27 

28 

29 

S7  8 

3 

3 

.3 

3 

3 

S8  = 

0 

0 

0 

0 

0 

S9  * 

0 

0 

0 

0 

0 

S10  ■ 

3 

3 

3 

3 

3 

Table  2:  Solutions  to  exanple  problem  (al) 
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B  * 

300 

295 

290 

285 

280 

min  MSRT  * 

.0838231 

.0859625 

.0878466 

. 0900560 

.0927486 

decision 

2 

2 

2 

2 

2 

S2 

3 

3 

3 

3 

3 

S3 

3 

4 

3 

3 

3 

S4 

4 

4 

4 

4 

4 

S5 

5 

4 

4 

4 

4 

S6 

21 

23 

23 

22 

21 

S7 

2 

2 

2 

2 

2 

S8 

0 

0 

0 

0 

0 

S9 

0 

0 

0 

0 

0 

S10 

3 

3 

3 

3 

3 

Table  3.  Solutions  to  example  problem  (bl ) 


B  - 

max  PA  * 

300 

.0403727 

29b 

.0387421 

B 

28b 

. 03b5446 

280 

.0341911 

decision  Sj 

3 

3 

3 

3 

3 

S2 

4 

4 

4 

4 

4 

S3 

5 

4 

S 

5 

b 

S4 

i> 

b 

b 

b 

b 

S5 

3 

3 

2 

2 

2 

S6 

26 

26 

26 

27 

26 

S7 

4 

4 

4 

4 

4 

S8 

0 

0 

0 

0 

0 

S9 

0 

0 

0 

0 

0 

S10 

2 

2 

3 

2 

2 

Table  4.  Solutions  to  example  problem  (cl) 

2.  The  performance  constrained  problems. 

The  related  performance  constrained  problems  were  also  solved  for  illus¬ 
tration.  For  example,  the  problem  (a2)  was  solved  using  the  data  from  Table  1 
with  the  restriction  that  SMA  >.  .732964  .  The  solution  to  that  problem  is  the 
same  as  the  solution  shown  in  the  last  column  of  Table  2  since  the  value 
.732964  is  the  (largest)  value  of  SMA  corresponding  to  a  budget  of  280. 
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III.  Modifications  to  Basic  Method 
A.  An  Alternative  Recursion 

An  alternative  and  more  efficient  approach  Is  available  for  the  budget 
constrained  problems. 

The  approach  Is  based  on  a  different  recursion  from  that  used  In 
Section  2  and  is  very  similar  to  an  approach  used  to  solve  the  "cargo  loading 
problem".  See  for  example  Dreyfus  [1]. 

The  cargo  loading  problem,  stated  as  a  maximization,  is: 

N 

max  l  v.d. 
j-1  J  J 

s.t.  I  c.d,  <  B 
j  J  J 

d-  ^  0  integer  . 

J  " 

Although  many  methods  are  available  for  solving  this  problem,  the  one  of 
interest  to  us  is  based  on  the  following  recursion 

f(b)  ■  max  {v.  +  f(b-c.)} 
je{l,...,N}  J  J 

where  f (b )  is  the  optimal  total  return  that  can  be  obtained  when  a  budget 
of  b  Is  available. 

To  Illustrate  this  method  consider  the  data  in  Table  5. 


1  1 

2 

3  “  N 

V1  1 

4 

6 

B  *  10 

ci  1 

3 

4 

Table  5.  Data  for  exanple  using  alternative  recursion 

The  computation  proceeds  with  Increasing  values  of  b  until  b  *  B  is 
reached.  The  process  can  be  viewed  as  shown  In  figure  3  where  a  template 
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representing  the  available  items  is  placed  over  the  budget  value  of  current 
Interest.  The  template  points  back  to  previously  determined  optimal  solution. 
Each  of  these  previous  solution  is  considered  for  updating  by  including  one 
more  item  of  the  type  inducted  by  the  template.  The  updated  solutions  are 
compared  and  the  best  selected  as  the  solution  for  the  current  value  of  b  . 
The  Illustration  shows  the  tenplate  at  the  budget  value  of  7.  The  optimal 
solutions  for  b  *  0,1,. ..,6  have  already  been  computed.  The  comparison  at 
B  *  7  is  among  the  solution  at  6  with  an  additional  item  1  for  a  total  return 
of  9,  the  solution  at  4  with  an  additional  item  of  type  2  for  a  total  return 
of  10,  and  the  solution  at  3  with  an  additional  item  of  type  3  for  a  return 
of  10.  Either  of  the  last  two  is  chosen  and  recorded  as  an  optimal  solution 


at  7. 


b  =  0  1  2 

d, (b)  *  0  1  2 
di(b)  =0  0  0 
d|(b)  *  0  0  0 


f  (B)  ■  0  1  2 


0  1111 

4  6  7  8  10 


1_J _ L 

8  9  10 


*an  alternative  optimal 
solution  at  b  =  6  is 

dj  “  0  dg  =  2  d^  -  0 


Figure  3:  Illustration  of  solution  method  using  the  alternative  recursion. 


A  very  simple  modification  of  this  procedure  can  be  used  to  solve  the 
problems  of  the  type  discussed  in  this  report.  We  consider 

n 


max  l  r.(d.) 
j-1  J  J 

n 

s.t.  I  c,  •  d,  <  B 

j-1  3  3 


d  >  0  Integer 

J 


where  the  return  functions  r.(d,)  are  concave.  The  procedure  below  is  not 

J  J 

guaranteed  to  give  the  optimal  solution  for  all  r.(d.)  but  is  guaranteed 

J  J 

if  the  r.(d.)  are  points  on  a  concave  function. 

J  J 

In  this  case  we  will  represent  the  return  functions  r.(d.)  as  the 

J  J 

sum  of  the  marginal  values  of  additional  items  of  type  j  . 

d  ' 

r,(d)  *  +  l  m,(i ) 

J  i=l  J 

Thus  r.(2)  »  m.(0)  +  m.(l)  +  m - (2)  .  These  marginal  values  form  a 

J  J  J 

sequence  with  the  properties  that 

and 

mj  ( i )  >  mj  lk )  i  <  k  . 

The  same  algorithm  as  before  was  applied  with  the  modification  that  the  value 
term  v,  ,  which  was  formerly  constant,  is  replaced  by  m.(i)  for  the 

J  J 

appropriate  value  of  i  .  The  program  which  implements  this  algorithm  was 
called  RECUR. 

It  should  be  noted  that  the  discussion  above  treats  the  constraint  as 
an  Inequality,  but  the  function  f(b)  In  this  section  is  computed  for  the 
constraint 

N 

l .  cjxj  *  b  • 

J»1  J  d 

For  this  reason  we  may  have  in  a  maximization  problem 

f(bj)  >  f(b2) 
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although  b^  <  b2  .  That  Is,  if  we  require  the  equality  to  be  met  exactly, 
it  is  not  necessarily  true  that  a  larger  budget  is  better.  The  program  prints 
the  values  of  f (b)  for  several  values  of  b  so  the  optimal  value  of  b  ^  B 
can  be  found  visually. 

B.  Minimum  Orders 

1)  RECMOD 

A  modification  was  made  to  the  program  RECUR  to  permit  the  user 
to  specify  minimum  packaging  quantities  of  each  item.  That  is,  item  i  is 
assumed  to  be  packaged  with  items  per  package.  The  provisioning  can 
select  only  whole  packages  of  each  item.  This  modification  resulted  in  the 
program  RECMOD  given  in  the  Appendix. 

2)  Exanple 

To  illustrate  the  RECMOD  program,  consider  the  example  problem 
(a2)  solved  previously.  The  optimal  solution  for  budget  B  *  300  is  repeated 
in  Table  2  for  the  case  in  which  all  q^  *  1  . 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Objective 

*1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

<i 

2 

3 

3 

4 

5 

21 

2 

0 

0 

3 

.0838231 

*1 

2 

1 

1 

5 

1 

1 

1 

3 

1 

1 

di 

2 

3 

3 

5 

4 

23 

2 

0 

0 

3 

.0847410 

Table  6.  Solution  to  Example  using  RECMOD. 


The  optimal  solution  is  also  shown  for  a  modified  problem  In  which  not  all  q^ 
are  equal  to  1  . 


16 


C.  Oiscussion 

The  modified  recursion  just  discussed  has  been  implemented  for  the  prob¬ 
lem  (a2)  which  is  to  minimize  MSRT  subject  to  a  budget  constraint.  The  method 
can  also  be  applied  to  the  other  budget  constrained  problems,  but  this  has  not 
yet  been  done.  All  that  is  required  is  to  modify  the  program  to  conpute  SMA 
or  PA  instead  of  MSRT  and  to  maximize  instead  of  minimize. 

There  Is  a  difficulty  in  extending  this  method  to  problems  in  which  the 
item  costs  are  arbitrary  values.  The  method  is  very  effective  when  the  costs 
are  all  integer  and  can  be  scaled  so  that  the  smallest  cost  is  1  .  If 
arbitrary  costs  are  allowed,  the  algorithm  can  become  ineffective  for  all 
except  small  values  of  B  .  Consider  for  exanple  the  costs  of  $l.uu,  $1.21, 
$1.27  for  three  items.  Let  the  budget  be  $2b.U0.  The  problem  could  be  solved 
by  scaling  the  costs  to  be  100,  120,  and  125  and  the  budget  to  be  2500,  but 
then  too  many  values  of  b  must  be  considered  when  many  of  them  are. not  pos¬ 
sible.  Alternatively  the  algorithm  can  step  to  the  “next  possible  value"  and 
will  consider  the  folowing  sequence  of  values 

b  *  100,  121,  127,  200,  221,  227,  242,  248,  254,  300,  321,  327 . 

As  the  process  continues,  depending  on  the  relative  values  of  the  costs,  the 
sequence  becomes  more  dense  and  may  eventually  include  all  possible  values  of 
b  .  This  Is  Ineffective  and  cumbersome  for  large  values  of  B  . 

It  is  also  not  possible  to  apply  the  modified  recursion  for  the  perform¬ 
ance  constrained  problems.  On  the  other  hand,  the  budget  constrained  problems 
are  solved  very  rapidly  and  the  relationship  between  cost  and  performance  can 
easily  be  determined  from  solving  the  budget  constrained  problem  using  the 
first  recursion.  In  fact.  If  MSRT  Is  minimized  for  a  budget  of  B  ,  the 
solution  Is  actually  obtained  for  all  values  of  b  up  to  and  Including  B  . 
Those  results  reveal  the  relationship  between  performance  and  budget. 


APPENDIX- PROGRAM-  LISTINGS 


FILE:  OPFILE  EXEC  A1  NAVAL  POST 3 R AD J ATE  SCHOOL 


FILEDEF  07 
FILEOEF  08 
FILEDEF  06 
FILEOEF  02 


0 ISK 
DISK 
TERM 
DISK 


(PERM 

RUSSIO 

(LAECL 

(PERM 


OATA 
133  R 


FM  FB  PE  *  M 


FILE:  FILE 


FTO7F00 1  A1  NAVAL  POSTGRADUATE  SCHOOL 


10 

NEW 

1 

10  0. 

3  00. 

300. 

300. 

10 

CLD 

2  95. 

295. 

10 

OLD 

290. 

2  9 0. 

10 

OLD 

285. 

2  85. 

10 

OLD 

280. 

280. 

10 

CLD 

2  75. 

2  75. 

MIN  SUM  2 


FILE:  RUSSIO  DATA 


Ai  NAVAL  POSTGRADUATE  SCHOOL 


2  0 


1: 

3. 

il¬ 

ls 

2. 


li 

10. 

2?: 

‘if: 

10. 


FOR  RECUR 
CURRENT  L 


30|(08,JON7050ST 


THE 

INE. 


15MA  ^MSRT^3AVA|L  - 

LAST  TWO  LINE  MUST  BE  PLACED 
IT  IS  FOR  DP4 


5 


.0137 

.0822 

}ih 

3082 
.1370 

DATA  ACE8DsLAM»C0ST  »  ES 
FIRST  INSTEAD  OF  THE 


»T,M7? 
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FILS:  0P4 
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OAT 4  DESCRIPTION  FOR  CMO  NJNBER 


columns'1 

justify 

V ARI BLE 
NAME 

MEANING 

l  -  5 

COL  5 

NMAX 

THE  NUMBER  OF  STAGES  IN  THE  PROBLEM 

6  -16"’ 

"col’IS" 

"  N  TAP  £  "  " 

LOGICAL  TAPE  NO  OF  THE  MAIN  TAPE 

IF  THIS  IS  LEFT  BLANK,  THE  COMPUTER 
WILL  USE  STANDARD  SCRATCH  TAPE  3 

II  -  15 

col  is 

SOLVE 

'-'  OLD  IF  THE  OPTIMAL  DECISION 

FUNCTIONS  HAVE  ALREADY  BEEN 
CALCULATED  AND  ARE  LOADED 

ON  LOGICAL  TAPE  NT  APE 

•  NEW  IF  THE  OPTIMAL  DECISION 

FUNCTIONS  MUST  BE 

CALCULATED  AND  STORED  ON 
LOGICAL  TAPE  MTAPE  BEFORE 
SOLVING  THE  PROBLEM 

DATA  DESCRIPTION  FOR  THE  NEXT  GROUP  OF  CAROS 

IF  SOLVE  «  NEW  ON  CARO  NUMBER  1.  TH=N  EACH  OF  T HE  NMAX  STAGES 
MUST  BE  DECRI BED  BY  THE  CARDS  DISCUSSED  BELOW.  YOU  CAN  USE  1  CARD 
PER  STAGE  OR  YOU  CAN  MAKE  1  CARO  DESCRIBE  MANY  ADJACENT  STAGES  IF 
THEY  ARE  SIMILAR.  THE  STAGES  MUST  BE  DESCRIBED  IN  NUMERICAL  ORDER 
STARTING  WITH  STAGE  1  II. E,  STAGE  1.  THEN  2.  ....  THEN  NMAX) 

OMIT  THIS  PACK  OF  CARDS  WHICH  DESCRIBE  THE  STAGES  IF  SOLVE  »  OLD 


COLUMNS 

£~"o~ 


1 1_-  20 
21  -’lo 

31_-  40_ 
41  -  46~ 


47  -  52 


JUSTIFY 

col's’*" 

col’Io” 


ANY 

•**««**« 

ANY 

ANY 

COL  46 


COL  52 


VARIBLE 

NAME 

NSTAGE 

’noTtto”" 


XLOW 

"xhigh 

oIlx’ 

XMOoi 


XSTAGE 


MEANING 


LOWEST  NUMBERED  STAGE  FOR  WHICH 
THIS  CARD  APPLIES 


HIGHEST  NUMBERED  STAGE  FOR  WHICH 
THIS  CARD  APPLIES.  IF  THESE  COLUMNS 
ARE  LEFT  BLANK,  THEN  NOITTO  WILL  BE 
TA<EN  AS  EQUAL  TO  NSTAGE__ 

lowIst’valuI  of *xn' for' stag!  nstagI’ 
highest’valuI  of'xn'for  stage  nstagI 
’incremInt  In'xn'for'stage  nstagI 

’*  "min’  IF  STAGE  NSTAgI  Is’tO  BE 
MINIMIZED 

’«  "max’  if’stage’nstage  li’To'ie  ’ 

^MAXIMIZED 

’«  "sum’  if’the’com  position* 

OPERATOR  BETWEEN  STAGES 
N STAGE  AND  NSTAGE- 1  IS 
ADOITION 

«  ’mult’  i  f’thI’ com posit  ion’ 

OPERATOR  BETWEEN  STAGES 
NSTAGE  AND  NSTAGE- I  IS 
MULTIPLICATION 


■HI  NMAX 


IF  THE  COMPOSITION 
OPERATOR  BETWEEN  STAGES 
NSTAGE  AND  NSTAGE- 1  IS  TO 
MINIMIZE  THE  MAXIMUM 
INDIVIDUAL  STAGE  RETURN 
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53  -  55 


COL  55 


MODES 


■NAXMIN  IF  THE  COMPOSITION 

OPERATOR  BETWEEN  STAGES 
NSTAGE  AMO  NSTAGE-1  IS  TO 
MAXIMIZE  THE  MINIMUM 
INDIVIDUAL  STAGE  RETURN 

■  this  varIable'is  not  used" 

WHEN  NSTAGE*!,  THEREFORE 
COLUMNS  47  TO  52  ARE 
IGNORED  WHEN  NSTAG E*i  AND 
MAY  BE  LEFT  BLANK 
HOWEVER  IF  NOITTO  IS 
GREATER  THAN  1,  THEN  XSTAGE 
MUST  BE  SPECIFIED  IN 
ACCORDANCE  WITH  THE 
COMPOSITOM  OPERATOR  FOR  THE 
REST  OF  STAGES  DESCRIBED  ON 
THIS  CARO 

~~I  I  F~ThI”fI BONACc7”siARCH  ll" 
TO  BE  USED 

IF" THE  FIBONACCI  sIaRCH  IS 
NOT  TO  BE  USED 


DATA  DESCRIPTION  OF  LAST  CARO  FOR  PROBLEM 

THE  OPTIMUM  DECISIONS  ARE  TO  BE  PRINTED  OUT  FOR  A  PROCESS  WITH 
NMAX  STAGES*  FURTHER,  FN(XN)  IS  TO  BE  OPTIMIZED  WITH  RESPECT  TO  XN 
WHEN  THE  PROGRAM  COMPLETES  THE  REQUESTED  PRINTOUT  IT  WILL  RETURN 
TO  THE  BEGINNING  AND  START  A  NEW  PROBLEM  IF  THERE  IS  DATA,  IF  NO 
DATA,  IT  WILL  EXIT 


COLUMNS 


1-15 


16  -  30 


JUSTIFY 


VARIffLE  •  MEANING 
NAME 


•  AT  STAGE  NMAX,  OPTIMIZE  FNtXN)  FCR 
.  XN  NOT  LESS  THAN  XM1 

AT  STAG  I  NMAX,"0PTIm7zI”fN(XN)"fCr" 

•  XM  NOT  GREATER  THAN  XN 2 


DESCRIPTION  OF  THE  SUBROUTINE  TRANFM 

THE  MAIN  PROGRAM  TRANSMITS  XN.DN,  AND  NUMBER  (THE  SEQUENCE  NUMBER 
OF  THE  PROBLEM!  TO  THE  SUBROUTINE  ANQ  THE  SUBROUTINE  CALCULATES 
YN  (THE  OUTPJT  OF  THE  STAGE  I.E.  YN*XN-1» 


DESCRIPTION  3F  THE  SUBROUTINE  STGRET 

THE  MAIN  PROGRAM  TRANSMITS  XN,  DN,  YN,  AND  NUMBER  TO  THE  SUBROUTINE, 
ANO  THE  SUBROUTINE  CALCULATES  RN  (THE  STAGE  RETURN! 


DESCRIPTION  OF  SUBROUTINE  STORE 

THE  MAIN  PROGRAM  CALLS  THE  SUBROUTINE  RIGHT  AFTER  READING  THE 
FIRST  DATA  CARO  AND  OEFINING  NUMBER.  THE  PURPOSE  OF  THIS 
SUBROUTINE  IS  TO  ALLOW  THE  STORING  OF  CONSTANTS  IN  COMMON  STORAGE 
FOR  THE  POSSIBLE  USE  OF  SUBROUTINES  STGRET  ANO  TRANFM 

DATA  WHICH  WILL  BE  RE AO  FROM  DATA  CARDS  AT  OBJECT  TIME  BY 

THE  SUBROUTINE  STORE  MUST  BE  INSERTED  BETWEEN  THE  FIRST  DATA  CARO 

ANO  THE  GROUP  OF  CARDS  WHICH  DESCRIBE  THE  INDIVIDUAL  STAGES 
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DESCRIPTION  3F  SUBROUTINE  OLIMIT  _ _ _  „  ,  „  _ . 

THE  MAIN  PROGRAM  TRAN  SMI  TS  N,XN»  AND  NUMBER  TO  THE  SUBROUT  INE, 

ANO  THE  SUBROUTINE  CALCULATES  DLOW  »  DHIGH,  AND  DELD 

WHERE  OLOW  *  LOWEST  VALUE  OF  ON  FOR  THAT  PARTICULAR  VALUE  OF  N 


ne  auBKjuiine  uuulaics  yu»»  unjum  ukuj 

■  OLOW  *  LOWEST  VALUE  OF  ON  FOR  THAT  PARTICULAR  VALUE  OF  N  ANO 
XN 

DHIGH  «  HIGHEST  VALUE  OF  DN  FOR  THAT  PARTICULAR  VALUE  OF  N 
AND  XN 

DELD  -  INCREMENT  IN  ON  FOR  THAT  PARTICULAR  VALUE  OF  N  ANC  XN 
IF  ANY  OF  THESE  VAR  ISLES  DO  NOT  CHANGE.  THEY  CAN  BE  SET  IN 
SUBROUTINE  STORE  AND  NOT  MENTIONED  IN  SUBROUTINE  DLIMIT 


IF  ANY  SUBROUTINE  SETS  NUFF  TO  -l  THEN  THE  PROGRAM  IGNORES  THE  _ 
VALUE  OF  ON  BEING  PROCESSED  AT  THAT  MOMENT ,  AND  IGNORES  ALL  LARGER 
VALUES  OF  DN  FOR  THAT  PARTICULAR  VALUE  OF  XN  AT  THAT  PARTICULAR 
STAGE — THIS  FEATURE  IS  DISABLED  4  HEN  USING  THE  FIBONACCI  SEARCH 


THIS  PROGRAM  USES  TAPES  3  AND  4  FOR  SCRATCH 
NSS2  ZERO  FOR  ENTIRE  PROBLEM  SOLVING  AND  OUTPUT 
NSS2  NON-ZERO  FOR  OUTPUT  FROM  PREVIOUS  PROBLEM 
DIMENSION  FNC 10001)  ,FNM1  (  10301 1 ♦  ONOFXNd 0001 ) 

REAL  *  8  HOLM  IN/*  MIN*/.  HOLOLO/  *  OLD*/ 

REAL  *  8  XMODE,  X STAGE.  SOLVE 
REAL  *  8  TAB)  4) 

COMMON  /C0N/TA6 

COMMON  XN.N.ON.  YN  .RN.  SLOW  .XHIGH.  DELX.  DLOW.DHIGH,  DELD,  NUMBER 
COMMON  YLOW. THIGH .DEL Y, NUFF, FN.JT OP, N NFL AG, NMAX.XLAM 
COMMON  A.B.C.O.E.SUMELT.  K00E1.K0DE2 
EQUIVALENCE  )FN)  1  )  *  FNMll  1 1  .ONQFXN  (III 
WRITE)6, 10001  HOLMIN.  HOLOLO,  (TABU),  1*1,4) 

1000  FORM AT (IX,  &A6I 

CALLPiEARCH(0.0,0.0,-l,0,0,0,DN3EST,BEST) 

NPAGE*0 

NUMBER-0 

240  REA0(  I NTAPE,  1 00  ,  END-2  000  )  NMAX  ,MT  APE,  SOLVE 
100  FORMAT ( 21 5,A5) 

MTAP  *  MTAPE 

NDITTO-O 

NCALC*4 

I F(M TAPE-4)  405,501,405 
501  NCAL*  * 


LC*4 
M TAPE-4) 

V8*  I72 

MTAPE)  3; 


GO  TO  172 

405  I  F(MTAPE)  322,322  ,172 
322  MTAPE* 3 
172  NTAP  E*MTAPE 
REWIND  MTAPE 
NUMBER-NUMBER+1 

CALL  STORE  WAS  MOVED  FROM  HERE 
I  FIS  OL  VE-HOLQLD )  301,302,301 
302  NSS2*1 

WRItI ( NOUTPE U52 )  NPAGE 
152  FORM  AT (1H1100X4HPAGEI 4) 

WRITE( NuUTPE,303)  MTAPE 


UirEi^v^,  nmrc 

58H_  f  HE  _  I  ABLE  S  QF  FN(XN)  _  AND  ONCXN^WWCH 


303  FORMAT  (58H  THE  TABLES  I 
1 APEI 5, 35H  WILL  BE  USED 
READ  (MTAPE)  NMAX 
REWIND  MTAPE 
GO  TO  304 
301  NSS2-0 

NpkbE-NpSlE+1 

WRITE) NQUTPE, 152)  NPAGE 
WRIT  Et  NQUTPE, 305 )  MTAPE 
305  F0RMAT(89H  THE  COMPUTER 


SOLVE  THIS 


LOGICAL 


TO  CALCULATE  TABLES  OF  FN(XN) 


DN(  XN 
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t it  AND  STORE  THEM  un  luu 
REWIND  NCALC 
500  WRITE! PTA PEI  NMAX 
304  WRITE ( NOUTPE,  3061  NMAX 
306  FORM  AT ( 16H  THE  PROBLEM  H 
314  IFINSS2)  197,109,197 
109  LINES  *  6 

WRITE£NQUTPE»2006 )  NMAX, 
2006  FORMAT!///  •  DATA  CARD 
121  DO  123  N-l.NMAX 

241 


LOGICAL  TAPE  I  51 


HASI5.7H  STAGESI 


MTAP,  SOLVE 
» ,  217,  A7) 


Mhhi 


241  READ (I NTA PE, 124  I  NST AGE, NDITTO, XL OW,X HIGH, DELX,XMODE, XSTAGE, MODES 

124  FORMAT (21  5,3c  10* 6  ,2A6 , 13  I  _ 

WRITE (NOUTPE,  2007 )  NSTAGE ,  N0ITT3,  XLOW,  XHIGH,  OELX,  XMODE, 

1  XSTAGE,  MODES 

2007  FORMAT!///  •  DATA  CARD  *,  215,  2X,  3F10.6,  2A6,  131 
I F(L IN ES—41 1  30,30,31 
31  NPAGE-NPAGE+1 

WRITECNOUTPE, 1521  NPAGE 
LlNES-1 

30  LINES  «  LINES  +  9 

I  TOP*! XHIGH-XLQW I /DcLX+1,  001 
IF(MODES-l)  600,601,600 

600  MODE  S--1 

601  1F!NSTAGE-N)  125.245,125 

125  WRIT E( NOUTPE ,  127 ) 

127  FORM AT ! 23H1DATA  CARD  OUT  OF  ORDER) 

GO  TO  2003 

245  I  F !  N  ST  AGE-ND I TTO  I  246,126,126  _ 

246  WRITE! NOUTPE,  247)  N.NSTAGE ,NDI  TTD 

247  FORM  AT ! ///46H  THE  DESCRIPTION  WHICH  APPEARS  BELOW  FOR  STAGE,I5,17H 
1  APPLIES  TO  STAGE, 15, 11H  THRU  STAGE,I5,10H  INCLUSIVE) 

LINES-LINES+4 
X XMODE  «  XMODE  -  HOLM  IN 
WRITE! 6,1001)  XXMODE 
1001  FORM  AT  (IX,  A6  I 

126  IF! XMOOE-HOLM  IN  )  230,231,230 
231  XMODE* 1*0 


W RITE ! NOUTPE, 232 )  N 

232  FORMAT  ( //6H  STAGE  15,1  8H  IS  A  MINIMIZATION! 

GO  TO  520 
230  XMOD E*-l. 0 

WRITE! NOUTPE, 2341  N 

234  FORMAT  !//6H  STAGE  I5.18H  IS  A  MAXIMIZATION) 

520  DO  129  1-1,4 
IF(XSTAGE-TAB(II I  129,130,129 

129  CONTINUE 

521  NNFLAG-5 
GO  TO  502 


>0  N NFL  AG-1 

^02  IF(ITOP-MAX)  316,316,112 
316  WRITE  (NOUTPE,  110)  XLOW.XHIGH.DaX 

110  F  ORM  AT  ( 63H  THE  STATE  VARIBLE  IS  BEING  TREATED  AS  DISCRETE,  GOING  F 
1R0M  XN-E14.5, 7H  TO  XN-E14.5.13H  IN  STEPS  OF  E14.5I 
IF(M0DES)T318^3ia  ,321 

lio  F0^Ma1?41H  t5tJl  *  ENUMERATION  IS  USED  FOR  THIS  STAGE) 

GO  TO  320 


ENUMERATION  IS  USED 


THIS  STAGE) 


>«,i  WRITE! NOUTPE t  , 

>11  FORMAT (40H  FIBONACCI  SEARCH 
>20  GO  TO  (11,12, 13, 14,1321, NNFI 
11  WRITE! NOUTPE,  16  I  NMl.N 
16  FORMAT (40H  THE  COMPOSITION  i 
IIS  AOOITIONI 


511) 


320  GO  TO 


RCH  IS  JS ED  FOR  THIS  STAGE) 
NNFLAG 


OPERATOR  BETWEEN  STAGESI6,4H  ANDI6,12H 


GO  TO  132 


12  WRITE! NOUTPE,  17)  NMl.N 

I7lI?RMUL4t?LlJi!flSN)P0^ITI0N  0PEMT0R  BETwEEN 


STAGESI6, 4H  ANDI6.18H 


GO  TO 


13  WRITE  (NOUTPE,  18 1  NMl.N 

18  F0RMAT!40H  THE  COMPOSITION  OPERATOR  BETWEEN  STAGESI6.4H  AN0I6.34H 
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TRAiVFM  TO  INDICATE  AN  INFEASIBLE  XN-1 


LIS  TO  MAXIMIZE  THE  MINIMUM  RETURN  I 
GO  TO  132 

WRITE (NOUTPE, 19)  NMi.N 

FORMAT  (40H  THE  COMPOSITION  OPERATOR  BETWEEN  STASESI6.4H  ANDI6,34H 
LIS  TO  MINIMIZE  THE  MAXIMUM  RETURN  I 
00  136  I»1,IT0P 
BEST  »1.0E+35*XM0DE 
X-I-l 

XN*X  LOW+X*D£LX 
NUFF*1 
CALL  OLIMIT 
I FtMODES )  602,602.603 

CALL  SEARCH!  OLOW,  DHIGH,DELD,XMQ3E  ,  DNB  EST  ,  BEST  ) 

GO  TO  400 

KTOP«<  OHIGH-DLDW ) /DEL D+l *001 

00  137  J«1,KT0P 

X«J-1 

DN*D LOW+X*DELD 

JUMP  MAY  BE  SET  TO  1  IN  TRANFM  TO  INDICATE  AN  INFEASIBLE  XN-1 

JUMP»0 

CALL  TRANFM 

CALL  STGRET 

XM0DE*1  FOR  MIN,  -1  FOR  MAX 
IF<  JUMP.EQ.il  RN*999999*XMOOE 
I  F(NUFF)  400,401,401 
IF(N-1I  504,504,503 
QN»RN 
GO  TO  143 

K*( YN-YLOWI/DELY* 1.001 
IF(K)  212,212,209 

rwiTH  X  N*E1  5.8 ,8H  AND  0N*E15.8,  12H  XN-1 
L EQUALS E16.8, 21H  AND  IS  OUT  3F  UNITS///) 

GO  TO  2003 

I F(  K— J TOP  I  205,206,21 2 
RNMl-FNMl(JTDP) 

-GO  TO  207 
X-K-l 

X»YLOW+X*OELY 

RNM1«FNM1(K)*(FNM1(K«-1I-FNM1(KII  MYN-XI/OELY 
GO  TO  (21.22,23, 24,248), NNFLAG 

FORMaI(43HiE5mPo!Aqn  OPERATOR  NOT  OEFINEO  FOR  STAGE, 151 

GO  TO  2003 

ONsRN+RNMl 

Sn»rn*)nhi 

tww  NM1)  141,141,142 

QN«RN 

GO  TO  143 

QN-RNM1 

GO  TO  143 

IF(RN-RNMl)  145,145,146 

QN-RNM1 

GO  TO  143 

QN-RN 

1 F  U  QN-BESTI  *XNOOEI  144,137,137 

BEST-QN 

PN6EST-0N 

CONTINUE 

WRIT  E(NCALC)  BEST  «DNBEST 
CONTINUE 


E15.8,  12H  XN-1 


141,141,142 


145,145,146 


144,137,137 


5Ifi|A8cftbSict>ON'D 

reaoTncalci  FN( 1 1 


tX^OW,XHIGH, 0ELX.IT3P ,XMQOE 
OUM,  ONOFXNi  I  I 


NCALC 


(ONOFXNt  II)  ,II-1,IT0PI 
TOP 

FN( I )  ,DUM 


23 


FILE:  DP4 


FORTRAN  A1  NAVAL  POSTGRADUATE  SCHOOL 


YLOW*XLOW 
YHIGH*XHIGH 
DELY-DELX 
JTOP-ITOP 
123  CONTINUE 

WRITE(NTAPE)  (FN(II),II» 
161  END  FILE  MTAPE 
197  REWlNO  MTAPE 

READ (MTAPE)  NMAX 

N-NMAX 

J*N-1 

IF(J)  407,407,201 


ITOP) 


irwi  tufttui 
00  522  1*1. J 
READ (MTAPE) 


522  READ (MTAPE) 

407  READ  (MTAPE)  XLOW  ,  XHIGH.OELX,  ITOP ,  XMOD E 
READ (MTAPE) 

R EAO (MTAPE)  (  FN (  I  ) ,  1*1 ,  ITOP) 

BACKSPACE  MTAPE 
BACKSPACE  MTAPE 
READ  (INTAPE, 162)  XN1.XN2 
162  FORMAT (2E 15. 8) 

NPAGE*NPAGE+1 
WRITE ( NOUTPE ,1521  NPAGE 
WRITE ( NOUTPE, 2008 )  XN1,  XN2 
2008  FORMAT!///  •  OATA  CARO  ',  2F20.3  III 
IF(XMOOE)  323,323,324 

3§t  F 0^514 (29HPfhll2?^0BLEM  IS  TO  MINIMIZE 
GO  TO  327 

323  WRITE! NOUTPE ,326)  N 

326  FORM  ATI  2  9H  THE  PROBLEM  IS  TO  MAXIMIZE 

327  WRIT  EC NOUTPE,  328 )  XN1.XN2 

328  FORMAT  (42H  XN  IS  TO  BE  CHQOSEN  OPTIM 


MINIMIZE 


STAGE  PROCESS) 


328  FORMAT  (42H  XN  IS  TO  BE  CHQOSEN  OPTIM ALL 
ID  XN-1PE14.5) 

171  IF(XLOW-XNl)  175.175,  176 

176  WRITE!  NOUTPE » 177  )  XLOW.XHIGH 

177  FORM  AT  (47H0THE  PROGRAM  ONLY  HAS  INFORMATION  ON  XN 
1 ANOE 16*8) 

GO  TO  240 

175  IF(XN2-XHIGHI  181,181,176 

181  IFIXN1-XN2)  178,178,176 

178  I  XI*  (  XN1-XL0W)  /DELX+1. 001 
1X2*  (XN2-XLQW  )/DELX*1.00  1 
B  E ST  *1 .OE+35*  XMOD  E 
00  182  J*IX1»  1X2 

IF((  FN(J)-BEST)*XMOOE)  183,182,182 

183  BEST *FN(  J  ) 

JSAVE-J 

182  CONTINUE 
X-JSAVE-l 

XN«XLCW+X*QELX  _ 

184  FORMATU^H  Opt? MALXxft|l?E14.  5.  WN  OPTIMA 
READ  (MTAPE)  (  ONOFXNC I  I ),  1 1*1 ,  ITOP  ) 

186  ggONOFXNJJSAVE) 

204  WRiVe(N0UTPE*188 ) 

188  FORMAT! 7 HO  N1  8X2HXN1  8X2HDN16X4HXN-1) 


TO  MAXIMIZE  AI 6, 14H  STAGE  PROCESS) 

HQOSEN  OPTIMALLY  BETWEEN  XN* 1PE14.5, 8H  AN 


BETWEENE16.  8  ,4H 


i?4lMALXxft2lJIl4.5,l6H 

(  DNQFXNCI  I),  11*1,  ITOP 


OPTIMAL  RETURN*  1PE14. 5) 


203  (NOUTPE, 189)  N,XN,ON,YN 

189  FORMAT!  17, 1P3E20.  5) 

NN-NN+1 

Haftm  1 

IF(N-l)  244,193,193 
244  GO  TO  240 


00  406  1*1,4  _ 

SISocStaS I^XlSw  ,  XHIGH,0E LX,  HOP.  X 
READ (MTAPE)  (  ONOFXNC I II, 1 1*1, ITOP I 
L*(  YN-XLOW  I/O  El  X+ 1.001 
IF(L)  214,214,215 


XMOOE 
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214  L«1 

219  NP1-N+1 

WRIT E (NOUTPE » 21 0)  NPi ,  XN »  ON,  YN 
NN-NN+7 

215  IFtL-ITOP)  218,216,217 

216  0N«0N0FXN(LI 
XN*Y  N 

GO  TO  196 

217  L»IT0P 

GO  TO  219 

218  XN«YN 
194  X»L-1 

X»XL0W*X*0ELX 

ON-ONOFXN <L  JM 0N0 FXNI  L*1 1  -0N0F XN ( L ) ) *  ( XN-X l/DELX 
196  CALL  TRANFM 

I  F(  NN-501  203,203,220 

220  N  PAG E*  NP AGE ♦ 1 

WRIT E( NOUTPE,  152 )  NPAGE 
GO  TO  204 


112  WRIT EiNOUTPE, 1131  MAX.ITOP 

113  FORMAT (72H  THIS  PROGRAM  LIMITS  THE  NUMBER  OF  DISCRETE  STEPS  OF  THE 
IS  TATE  VAR  ISLE  T0I6,26H  AND  THIS  PROBLEM  REQUIRES  16 1 

2002  FORMAT  ( *1'  ,*  ERROR  HALT') 

2003  WRITE! 6, 2002) 

STOP 

2000  WRITEI6.2001) 

2001  F ORM AT I • 1  * , •  END  OF  DATA  FILE*) 

STOP 

END 

CSEARCH  DISCRETE  FIBONACCI  SEARCH  SUBROUTINE 

SUBROUTINE  SEARCH  IAA,  BB,  OELYY,  XXNODE,  YYBEST,BEST ) 

DIMENSION  FI150) 

OPTIMIZE  WITH  RESPECT  TO  Y  BETWEEN  AA  ANO  BB  IN  STEPS  OF  DELYY 
STORE  OPTIMUM  Y  IN  YY BEST  AND  TN5  OPTIMUM  VALUE  OF  THE  OBJECTIVE 
FUNCTION  IN  BEST 
XXMQDE»-1  FOR  MAXIMIZE 
XXMODE-1  FOR  MINIMIZE 

SUBROUTINE  MUST  BE  INITIALIZED  BY  CALLING  IT  WITH  DELYY*-1.0  AT 
LEAST  ONCE  BEFORE  IT  IS  USEO  FOR  A  SEARCH 
DELY-OELYY 
A«AA 

b«bb 

X  MOD  E* XXMQOE 
IFJDELVI  100,100,101 

100 


P?ll«i(l-1)2F(I-2)*1.  0 

IF(F(I)-i.0E+35l  1,2,2 
2  I  I«I 

1  ^SSjfiNue 

RETURN 
101  YHI-B 
YLO-A 

FNO- I YHI-YLO) /OELY+1. 0 
.  ??<FNOrMN))  6,6,5 


FORM A^ tSBHIERROR,  TOO  MANY  POINTS  TO  BE  SEARCHES! 
FORMAT!  ERROR  HALT') 


288  sw:^)6**0* 

STOP 

6  IF<N-2)  42,41,40 
40  Y l*r ( N-2 ) *DEL Y*  YL  0 
CALL  FUNCTN?Y1,GY1) 
Y2*FIN-1I*DELY*YL0 
CALL  FUNCTNl Y2,GY 2) 
16  N*N“1 
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IF(N-1»  102.102.103 

103  IF((  GY2-GY1 )*XMuDEl  19.19,20 

19  YLQ»  Yl+DELY 
Y1*Y2 
GY1«6Y2 

Y2*F(N-U*0ELY+YL0 
IRY2-8I  104,104,  105 
105  GY2«1.0E435*XM0DE 
GO  TO  16 

104  CALL  FUNCTNt Y2, GY  2) 

GO  TO  16 

20  YHI* Y2-DELY 
Y2«Y1 
GY2»GY1 

I  FIN-2 »  34,34,35 

34  FNM2=0.0 
GO  TO  36 

35  FNM2«F<N-2I 

36  Y 1*FNM2*QELY* YLO 
CALL  FUNCTNI Y 1 ,  GY  1 1 
GO  TO  16 

102  IF(<GY2-GYll*XM0DE>  110,111,111 

110  YYBEST-V2 
BEST«GY2 
RETURN 

111  .YYBEST*Y1 
BEST*GY1 
RETURN 

41  Y1*A 

CALL  FUNCTNm.GYll 
Y2*B 

CALL  FUNCTN(Y2,GY2) 

GO  TO  102 

42  Y 1*A 


CALL  FUNCTNC Yl.GYll 
GO  TO  111 


CFUNCTN 

SUBROUTINE  FJNCTNIY,SY» 

DIMENSION  FNI 10001)  *FNM1  (  1Q0Q1),  ONOFXNdOOOl 1 
REAL  ♦  8  TAB(4» 

COMMON  /CON/TAB 

COMMON  XN.N.DN,  YN  ,RN,  SLOW.XHIGH.OELX,  DLOW.DHI  GH, 
COMMON  YLOW.YHIGH.OELY.NUFF.FN.JTOP.NNFLAG.NMAX, 
COMMON  A.B.C.D.E,  SUMELT.  KQDfe 
EQUIVALENCE  i  FN(  1 1  ,FNM1(  1 1  ,0N0FXN(  111 


DELD, NUMBER 
XLAM 


ON«Y 

N0UTPE»6 

JUMP  MAY  BE  SET  TO  1  IN  TRAN 

JUMP«0 

CALL  TRANFM 

CALL  STGRET 

XM0DE*1  FOR  MIN,  -l  FOR  MAX 
IF( J  UMP,£Q, 1 j RM»99999 9*XM0DE 
I F(N~l I  146,146,401 
K*<YN-YL0WI/DELY«-1.001 
IFJKI  212,212.209 
WRITE  (NOUTPE*  210 I _N,XN,0N , YN 


IN  TRANFM  TO  INDICATE  AN  INFEASIBLE 


5^8, 8H  AND  DN«E15.8, 12H  XN-1 


r  Unfl  All*,*  CnKUn  HAkl’l 

>003  WRITE (NOUTPE, 2002 1 

209  IFiK-JTOPI  205,206,212 
iob 

205  X*K-1 

x*ylow+x*qely 

2ot  sn3f?si5i:iF3^L,^yiAFr<K",<w-<'/OEU 

n  mrw1 
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22  QN«RN*RNM1 

23  ?F(RN-InH1I 

141  137 

142  MIPilT,, 

i2!  ths-*ru 

GO  TO  137 
146  QN-RN 
137  GY-QN 
RETURN 
END 


141*141*142 

145*145*146 
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CSTORE 


/CON/TAB 


808 


110 

101 

202 

303 


COMMON 
COMMON 
COMMON 
COMMON  a 9 o • 
EQUIVALENCE 
I NDATA»8 

R?XJPf  IN0ATA«808I 
FORMAT  (2131 


(  10001), )NOFXN(10001I  ,A(101»  t B(  10  1 1.C (  101 1 


DHIGH#D ELD » NUMBER 
NMAX»XLAM 


KOOElt  K00E2 


U5881 


ifcffcfl 

IFIK0DE1.EQ.2) 

IF(K00E1.EQ.3I 


WRITEJ  NOUTPj 


ORM  AT 
FORMAT 
FORM  AT 
FORMAT 
I  F(KODE 
IFIKOO 
I  FtKOD 


WRITECNOUTPE.IOiI 
WRITEi NOUTPE. 202) 

MR  ITE(  NOUTPE 1 3031 
(3 OH  THE  OBJECTIVE  IS  COST  ) 

(30H  THE  OBJECTIVE  IS  SNA  I 

(3 OH  THE  OBJECTIVE  IS  MS*T  I 

(36H  THE  OBJECTIVE  IS  AVAILABILITY ) 

2.EQ.0I  mRITEI  NOUTPE, 900) 

2.EQ.1I  WR|T§(NOUTPi.?OII 
2 • EQ* 2 


WRITE! NOUTPE. 902) 


903 


999 

333 

30 

765 


FORMAT 
00  333 


AVAILABILITY 


WRIT 

CONT 


1 3 OH  CONSTRAINT 
ItNMAX 
ATA.99 
(5F10.4) 

.999)  A(K),B(K),CUl  >D(K)  ,E(K) 


A(K),C(K),E(K)  >B(K)  >D(K) 


!nu?utpe 

tTI»!,NMAX 


So'li1 

FORM  AT  (F10«3) 


RETURN 
END 
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CSTGRET 

SUBROUTINE  STGRET 

01 ME NS  ION  FN(10001I.FNM1(  lOOOllONOFXNdOOOll »A( 1  011 , B( 10 1 l»C( 101) 
DIMENSION  0(101 1  1 1(101 1 
REAL  *  8  T  AB  (  4 1 


COMMON  /CON/TAB 
COMMON  XN»N»3N*YN«RNf  SLOW  tXHIGH*  OELXy  OLOM.OHIGH,  OELO.  NUMBER 
COMMON  YLOW1YHlGHtOELVtNUFFiFNtjrOP*VNFLAGfNMAXfXLAM 
COMMON  A.BtC.DtE.5UMELT.K0Dcl*K0DE2 
EQUIVALENCE  (  FN(  1  I ,  FNMll  1 1  tONOFXM (  111 
IFUODE1.EQ.O  IGO  TO  77 
ION* INT(DN) 


A 8* A (N)*B(N) 

I  S*I  ON 
TERM*0* 

T  EMP*0« 

IF(IS.LT.O)  GO  TO  11 

TERM*EXP(-ABI 

TEMP-TERM 

I F( I S«  EQ.0I  30  TO  11 
00  10  1*1 »  IS 
TEMP-TEMP*A8/I 
IF(TERM.GE. .999991  GO  1 
IF(TEMP.LE.. 00001  I  GO  1 

10  TERM-TERM  +TEMP 

11  CDF-TERM 
20  CONTINUE 

IFUOOEl.NE.llGO  TO  40 


C  F  LLSMA??AB*flf-?EMP?i?ION-ABI*(l.-COFI)/AB 

sma* ( l •—Tempi  ♦( i on-ab i *( i.-cdfi/ab 

R  N*E ( N I  *AB*SM  A/ SU MELT 
RETURN 

40  TWUS*(  1,-CDFI  *(  AS*AB-2.*AB*I0N*I  0N*(  I  DN*1 1 1/  (  2.  *A  ( N)  I 
X*,TEMP*B(N)*(AB-IDNI/2  • 

IF(K0DE1.EQ.3|G0  TO  50 
RN-E (N I4TWUS/ SUMELT 


RETURN 


50  CONT 


AMSRT-TWUS/AB 

AMT8F*1./A(NI 

AMTTR-O(N) 

RN-( AMTBFI/ ( AMT BF+AMTTR*AMSRT I 

RETURN 

RN-C (N l*0N 

IF(K00E2.  EQ.l.AND.N.EQ.l.  AN0.YN.LT.XNIRN-999999 

RETURN 

END 
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COLIMIT 

SUBROUTINE  DLIMIT 

DIMENSION  FN(  10001)  *FNM1  (  10001 )»  ONOFXN(  10001  )«A(101),B(101)«C(  101) 
01  ME  NS  ION  0(101)  ,§(101) 

REAL  *  8  TABU) 

COMMON  /CON/TAB 

COMMON  XN,N,  ON,  YN  ,RN,  SLOW  ,XHIGH,OELX,  OLOW,OHIGH,D£LD,  NUMBER 
COMMON  YLOW, YHI  GH  ,OEL  Y,NUFF (FN*J  TOP,VNFLAG,NMAX,  XLAM 
COMMON  A«8,C,0«E, SUME LT, KGDel ,  K3DE2 
EQUIVALENCE  (  FN(  I  I.FNMli  1 )  ,DNOFXN (  l) ) 

IF(KODE2.NE.O)GO  TO  3  7 
DLOW-O. 

DELD-1. 

0 HIGH* AMI  N1  ( XN/C (  N)  ,50.) 

RETURN 

37  AB»A(N)*B(N) 

DO  555  10-1,30 

IS-IO-l 

TERM-EXP(-AB) 

TEMP-TERM 

IF(IS.EQ.O)  GO  TO  11 
00  10  I«1,IS 
TEMP»TEMP*AB/I 
IF(TERM.GE. . 99999  )  GO  TO  11 
C  IF(TEMP.LE.. 00001 )  GO  TO  11 

10  TERM-TERM  +TEMP 

11  CDF- TERM 

40  TWUS-( l.-CDFl  *(  AB*AB-2.*AB*I S* ISMIS*1))/(2.*A(N) ) 
X«-TEMP*B(N)*(AB-IS)/2. 

IF<KOOE2.EO.2.AN0.ETN)*TWUS/$UMELT.LE.XN)GO  TO  556 
C  CHANGED  FOLLOWING  TO  555  FROM  556 
IF(KODE2.EQ.2)GO  TO  555 


IF(KODE2.EQ.2)GO  TO  555 

IF(KODE2.EQ.3)GO  TO  554 

SSS-  ( l.-TEMP)  ♦(  IS-AB)  *(1.-CJF|/*B 

RR-E (N)*AB*SSS/SUMELT 

Y-XN-RR 

OLOW-O. 

THE  NUMBER  IN  THE  NEXT  LINE  HAS  A  BIG 
I F(N«  EQ«  1 IDLOW- 1000. 


EFF  ECT  . 


SHOJLD  NOT. 


I F  ( N 
IF(N 


.EQ.l.AND.Y.GT.O.  .AND. I3.LT.30 )G0 
•  EQ.l.AND.Y.LE.  0.  IDLOW-IS 


136  0HIGH-30.+0L0W 

I F(0 LOW.EQ.l 000.  I  DHIGH-1000. 


m*' 

554  AMSRT-TW 
AMTBF-1. 
AMTTR«p( 
AV-AMTBF 

iF(AV.GE 
ONTINUE 
S-30, 
LOW -IS 


TWUS/AB 

1./A1N) 


lAMTBF+AMTTR+AMSRT) 
(N)  GO  TO  55  6 


.  ILO-I. 

D  HIGH- 50. 

(KU"N 

BLOCK  OATA 
COMMON  /CON/TAB 
REAL  *  8  TAB( 4) 


MULT'  ,'MAXMIN'  ,  •  MINMAX'  / 


30 
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SUBROUTINE  TRANFM 

SlfflffllioN  5J{yj0^j|5jyi,10001,»5NOFXNl10001*  101> 

REAL  *8  TAB<4)  * 

COMMON  /CON/TAB 

g858!§a  ?£?6t3:?KiSS:g&fiSa^&”i?i^S^?I3«;3EkS*MOHBER 

EQU?  VALSVCl  I FN<  I  ^PN^ll  1 > ?Dn6fx3  <  I II 
I  F(  K0DE2.  NE*0  )60  TO  12 
YN*XN-C(N)*DN 


RETURN 

12  IDN*INT(DN) 

AB*A(N)*B(NI 

IS-IDN 

T£RM*0. 

TEMP*0. 

1F(I S«LT*0)  GO  TO  11 

TERM*EXP<-ABI 

TEMP -TERM 

1FUS.EQ.0I  GO  TO  11 
DO  10  1*1, IS 
TEMP-TEMP*AB/I 
IFITERM.GE..99999I  GO  TO  11 

i F(T Imp •!£•• ooooi )  go  to  11 

10  TERM*TERN  ♦TEMP 

11  CDF*  TERM 

40  TWUS*U«-CDF)  *t  AB*AB~2**AB*I0N*l  0N*(  I  DN+1 1  l/(2.*A(N>  I 
X+TEMP*8<N)*<AB-IDN)/2. 

I  F(KQDE2.  ED.  2)  GO  TO  50 


f-tiS?? ♦€ ?8n^AB ?*(  1«-C0F|F  AS 
YN*XN-E(N)*AB*SMA/SUMELT 


IFIYN.LT.O)  YN*0.0 


IFIN.NE.l.AND.YN.LT.O.OA.N.NE, 
IF(N*EQ,l*ANO.YN. GT.O )JUMP*1 
IFiJUMP.EQ.DYN-O, 

RETURN 

51  AMSRT*TWUS/AB 
AMTBF«1./A(N1 

amttr*o<ni 

AV*(  AMTBFI/<AMTBF+AMTTR*-AMSRT) 

YN-XN/AV 

AV  IS  AVAIL 

RETURN 

50  YN*XN-E(N)*T«US/SUMELT 
IF<YN.LE.0.)YN*0. 

RETURN 


(.LT.O.OA.N.N 
UGT.O  I  JUMP*! 


0  «0R«  N«N£*1» AN0*YN»6T*i« ) JUMP* l 


uo 
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DIMENSION  A(  20 )  «  B  (20t  ,CC  20 1 .  DC  23  I  •  E(20l,  MINCH  201 

DIMENSION  T(10,30),F<  1200  )  ,DX C  ID  ,  1203  )  ,TN(10 .  31 1 

INTEGER  BUDGET »  BB  »BBP  1  ,DX 

REAL*8  TIMEtCTIME 

I NDA TA*8 

N0UTPE*6 

READ  (INDATA, 8  88)  NMAX  ,  BUDGET 
888  F0RMAT(2I 5) 

READ  (  INOATA, 7TT )  (MINQ( I)«I*1«VMAX) 

7T7  F ORM AT ( 2015) 

DO  333  K»1 , NMAX 

READ  (  INDATA*  999 )  A(K)  ,C(K),E(K),B(K)  »  D(  K) 

999  FORMAT  (5F10.4) 

C  WRITE ( NCUTPE,  999 )  A(K  )  ,B  < K ) ,  C< K)  ,  D <K) , E< K ) 

333  CONTINUE 

TIME*CTIME(1) 

SUMELT-O. 

DO  30  1*1  ,NMAX 

30  SUME  LT*SUM£LT  •*■£(! l*A( 1 1*8  4  I) 


14 

11 


DO  10  1*1, NMAX 
OX(I  ,11*0 
DO  110  KP1*1, 31 
K*KP1— 1 
A B*A  ( I  ) *B( I ) 

TERM-EXP<-ABI 
TEMP*  TERM 

IF(K.EQ.O)  GO  TO  11 

DO  14  JI*1,K 

TEMP*TEMP*AB/JI 

IF4TERM.GE.. 999999)  GO  TO  11 

I F (TEMP* LEs* 00000 01 )  GO  TO  11 

TERM*TERM  +TEMP 

CDF*TERM 


40  7WUS*(1.-CDF>*(A8*AB-2.PAB4K>K*(K+1)I  /(2.*AII  )) 
X+TEM P*B(I)*(AB-K)/2. 

110  TN(I.KP1)-E(I)*TWUS/SUMELT 
DO  70  K*1 ,30 

70  T(I,KI*TNU,K+1)-TN(I  ,K» 

10  CONTINUE 


PQ.456 . 1  JKfl  j 


C  __  .... 

C  456  WRITE(6,555m  i 
C  555  FORMAT  (5F12. 7) 


KlfIJK),  K  1*1,  NMAX) 


90 

85 

80 


DO  80  1*1, NMAX 
LIM-30/MINQU  ) 

DO  85  K*1,LIM 
TMO0*0. 

MM*M INC(  I ) 

DO  90  J*1,MM 

TMOD»TMQD«-T(I  ,(K-1)*MINQ<  IM-JI 
T(I,KI«TMOD 

C?II-M^NQ(I)*C(II 

CONTINUE 

I BD-BUDGET+1 
DO  100  BBP1*1 ,1 BD 
BB*BBP1-1 
XMIN*0. 

ISTAR-0 
TRV*0* 

DO  210  1*1, NMAX 

iHItt1'*1-9' 

TRY-AMIN 
GO  TO  21 


GO 


210 

211 
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IMS!* 

210  MfflU NIN 

19  $  f  1*  8S  :il^|  J  i  IT  SS  i  :tl :  8 !  I8W,,# 

Aill!lIpil.DX(I.BiPl-CCISTARn 


19  ??  HsJUl§Sto.QR.gg-CjI 
ox(i!a8Pi)iox?i?BEPi-Eii 

“  IfeffiKff! 


iillls&mar.sn . . 

;nue  _ 


DX(I$TARtB8Pl  l*AHINOi  OX%  I5TAK*t»Fi"v.i  * 

100  CONTINUE 
KLOW*I BO-5 

??m|2+i51Hc?VS|?|Vigh 

°3456  5S??|H!li!i!5x(!:KBDIfI«l.NMAX» 

WRiTc  ( 6«l600l  iHINsI  I)  iNMAXi 
1000  FORNATjISlSl 

88  lUt-iSUm  I  ♦  « »  *0**  1  *  » +U 

WR?T l**6?2000 1  AMSRT* KB0M1 
2000  FORMAT <1X»*MSRT*  •,F10«8#*  BJ3GET- 

101 

?A?SteaO  I  GO  TO  10 

ift  Hi 


GO  TO  10 


HALLESETIMi 
R  iTURN _ 


10  l^EGH?fU?18!U26DO 

RETURN 
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